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We show that in magnetic nanostructures with a homogeneous magnetic order, the equilibrium 
spin current can be nonzero. For example, this is the case of a wide magnetic ring with the mag- 
netization along the ring axis. The physical reason of this effect is a variation of the orientation 
of anisotropy axis inducing a spin torque acting on the magnetic ions. The mechanism of the spin 
current generation is related to the quantum vacuum fluctuations in the magnetic system. 
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One of the key problems in modern magnetoelectron- 
ics is related to the creation and manipulation of the 
spin currentsi'^'^'^. It was shown recently that the spin 
current can be generated using the flux of magnons in 
a nonequilibrium state of the magnetic system^. This 
idea was further developed for the magnons in textured 
magnetic structures like magnetic rings in the inhomoge- 
neous external field^i^. In the latter case, the existence 
of spin currents is related to a gauge field for the motion 
of magnons^i^. The gauge field and, correspondingly, the 
Berry phase for the adiabatic motion of magnons appear 
naturally due to the nonvanishing chirality of the mag- 
netic ordering. In terms of the Berry phase, the magnon 
mechanism of the spin current excitation can be better 
understood for various types of magnetic structures^. 

In the case of magnetic systems with topological ex- 
citations (like vortices or skyrmions) or a system sub- 
ject to inhomogenous fields, there appears an equilibrium 
spin current associated with the transmission of angular 
momentum^. This current transfers the torque acting 
on the disoriented magnetic moments. The equilibrium 
spin current is usually much stronger than the one trans- 
ferred by magnons because it does not require any spin 
wave excitations. This mechanism of spin current gener- 
ation can be realized only in textured ferromagnets, and 
it disappears in case of a homogeneous magnetization. 

In real magnetic structures, one should also take into 
account a possible inhomogeneity with respect to the ori- 
entation of anisotropy axes. It induces a new type of a 
topological Berry phase, which does not necessarily van- 
ish when the usual geometric phase is exactly zero^. An 
example of such a system is a wide magnetic ring (thin- 
wall cylinder) with the magnetization along the cylinder 
axi o^'"*"" . In this system, the gauge field is related to the 
variation of anisotropy axis (perpendicular to the cylin- 
der surface). 

We show below that the equilibrium spin current does 
not vanish in the system with the homogeneous magne- 
tization but inhomogeneous anisotropy. This effect has 
purely quantum origin and is related to the vacuum fluc- 
tuations in the magnetic system. The quantum character 
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FIG. 1: Two spins S = \ with different orientation of 
anisotropy axes m x — y plane. 

of this effect makes it substantially different from the ear- 
lier studied mechanism of the equilibrium spin current in 
noncoUinear ferromagnets^jii, which can be understood 
within a model of the classical moment. 

The spin current appearing in a shaped mesoscopic 
structure (e.g., a nanoring) with the colhnear magnetic 
ordering, transfers the angular momentum. The resulting 
transferred torque acts on the anisotropy axes, tending 
to arrange them homogeneously. In other words, the spin 
current produces a mechanical stress in the anisotropy- 
noncoUinear system. 

The physics of this phenomenon can be better under- 
stood if we consider first a simple model with two quan- 
tum spins 5 = 1 in the magnetic field B directed along 
the axis z, and with the anisotropy axes oriented along 
rii and n2 in the x — y plane for spins S\ and 5*2, respec- 
tively (Fig. 1). We take ni along x and denote the angle 
between ni and Wi by 0. The Hamiltonian is 

n = -BSi, + ^ (m • Sif-BS2^ + ^ (n2 • S2)'-.gSrS2, 

where A > is the constant of anisotropy, and the last 
term corresponds to the exchange interaction with the 
coupling constant g. We assume g > 0, which corre- 
sponds to the ferromagnetic coupling. Obviously, in the 
ground state both spins are oriented along z. 

We can use the basis of eigenfunctions of Hamilto- 
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nian (1) with g = 0. For each of the noninteract- 
ing spins Si and 5*2, the cigcnfunctions can be easily 
found, and they differ from each other by the trans- 
formation exp(i0Jj), where is the generator of z- 
rotations. In the absence of interaction, g — 0, the low- 
est four states of Hamiltonian (1) correspond to energies 
(2£o, So + A/2, A), where the level Eq -|- A/2 is twofold 
degenerate. We denote these states |00) , |01) , |10) , 
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Here Eq = A/4 — (A^/16 + i?^) is the ground state 
energy of a single spin. 

The essential point is that each of two separated quan- 
tum spins in the ground state has the magnitude of the 
moment along the z axis smaller than unity, |(S'i2)q| — 
\{S2z)o\ = 2\eo\B/ {el + B^) < 1. This is the effect 
of the magnetic anisotropy. As a result, the quantum 
fluctuations give a different value of the fluctuating mo- 
ments along the axes x and y (principal axes for 5*1), 
(Si), - (Sfy)^ = {el - B^) I {el + B^). It leads to the 
dependence of interaction on a mutual orientation of the 
S\ and 5*2 anisotropy axes. 

The interaction Hint = ~.9Si • S2 gives nonzero matrix 
elements 
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Assuming i? 3> A, we can restrict ourselves by consider- 
ing only the four lowest energy states. Then the ground- 
state energy £'0 of the Hamiltonian (1) can be found by 
the diagonalization of matrix 4x4, and we obtain 
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where Eq — e^ + V/2. Due to the exchange interac- 
tion, the energy depends on the mutual orientation of 
the anisotropy axis. The torque acting on the spins is 
A^((/)) = —dEo/dcj). In the limit of g — > we find from 
(4) 



52 {el - BY sui2(t> 
'2{el + B^f{X/2~eo) 
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The torque vanishes for parallel or antiparallel configura- 
tion of anisotropy axes, and also if they are perpendicular 
to each other, (j> = 7r/2. 

Let us consider now the model with arbitrary spins 
S > 1 on a lattice, with a nearest-neighbor exchange 
interaction. We assume a one-site anisotropy with the 
anisotropy axis smoothly depending on position. The 
Hamiltonian has the following form 
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We assume that in the ground state the ferromagnet has 
a uniform magnetization along the axis z, and all vectors 
rii lie in the x — y plane. It corresponds to the easy plane 
ferromagnet with a nonuniform anisotropy and can be 
realized, for example, in a thin-wall magnetic cylinder 
(Fig. 2)8'io. 

Using a local frame with the x axis along in each 
lattice site, one transforms (1) into 
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where (jji is the rotation angle, and we use a ma- 
trix form of presentation for the spin operators, Sf = 

After the Fourier transformation, — 5^ g''^ ''' , 
and taking the limit of g ^ the Hamiltonian is 

^ = i E [^^q + 2^(q • (^-q^q - ^-q^q) 
q 

+ (5Vq + ^Vq)] + ^ E(^f )^ (8) 

i 

where A. = ^icjji is the gauge potential, a = gal, ao is the 
lattice constant, and we assume that the variation of ^ 
in space is small at the wavelength of magnons (adiabatic 
approximation )S. 

We use the Holstein-Primakoff representation of spin 
operators-12, = S - n„ S+ = V^S {I - nj2Sy^^ a„ 

= V2S (1 ~ ni/2S)^^^ , rii = a^a^, and restrict 
ourselves by the harmonic approximation in ai , . Here 
al , Qi are the boson creation and annihilation operators 
for the magnons. Then we find 

^ = f E { [(« + '^)' + ^1) ("q'^q + 
q 

+A2 (aqa_q -I- a+atq) } (9) 

where Ai = A (1 - 1/25) and A2 = A (1 - 1/25')^/^ 

In the following, wc make a change in Eq. (9) sub- 
stituting A2 ^ Ai = A (1 — 1/25"). This is related to 
the well-known renormalization of the constant A2 due 
to the magnon interactions^^. It reproduces correctly 
both S" — > 1/2 and 5 — > 00 limits, and leads to the 
gapless spectrum of magnons in correspondence with 
the Goldstone theorem^^. The Hamiltonian (9) with 
A2 — > Ai can be diagonalized using the Bogolyubov- 
Holstein-Primakoff transformation mcthod^^, and wc ob- 
tain finally 
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where 6^, 5q operators are related to a^, Oq by the 
Bogolyubov-Holstein-Primakoff transformation 

feq = aqflq + /?qa-q, fojj = «* flj^ -|- /3*a_q (11) 
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and Wq is the energy spectrum of magnons, 



UJq^ S (^Jq - J-q]+S 



Here we denoted 



Jq + J-cl][Jq + J-q + 2Al 
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J{0)^J{q)^a{q + Ay (13) 

and used the standard notation J(q) for the Fourier 
transform of exchange interaction. 

Using (12), (13) we can find the dependence of momen- 
tum q on A. It determines the Berry phase acquired by 
the spin wave, which moves in the gauge potential re- 
lated to the varying anisotropy. We find that in case of 
adiabatic motion (see details in Ref.^), the Berry phase 
along a contour C is 7b (C) = A(r) dr, where A. = 



A 



-1/2 



is the effective gauge potential 
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acting on magnons and — 2S [aq^ {aq^ -f Ai)] . For 
Ai 0, we get A = A but for a large anisotropy, 
Ai ^ the effective potential A ~ Aui^/SXi is 

small comparing to A. Thus, the Berry phase of magnons 
is determined by ^(r). It generalizes the result of Ref,^ 
for arbitrary quantum spins on the lattice. 

We consider now the contribution to the total energy of 
the vacuum fluctuations, Eq = X]q('^q/2)i corresponding 
to the second term in (10). The energy i?o depends on 
the gauge potential A resulting from the inhomogeneous 
anisotropy. To calculate the spin current density j" we 
add to A an additional fictitious potential A!{v)^ which 
we put zero after calculation. Thus, in the previous def- 
inition of Jq we substitute A A + A' . Note that 
A , like the original potential A, is associated with the 
spin transformation - rotation around the axis z, which 
corresponds to the spin current density with the spin 
polarization along z. 

We use the definition jf (r) = 7 [SEo/dAi{r)]^'^^ fol- 
lowing from the connection between the spin current con- 
servation and the invariance with respect to rotations in 
the spin spaced, where 7 is the gyromagnetic ratio. Then 
we obtain 
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where n is the volume. In the usual bulk or an infinite 
2D system, the first term vanishes due to the inversion 
symmetry q — > — q, and, taking the limit of ^ — > 0, we 
find 

2ag2 + Ai 
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Here the sum over momentum runs up to qmax — t^/o-o- 
We can estimate it for the 3D case as 
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FIG. 2: Thin magnetic cylinder with the homogeneous mag- 
netization and noncoUinear anisotropy. 
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where k = (Ai/a)^/^. 

As follows from (15), the spin current flows along the 
gradient of the anisotropy axis variation. In terms of the 
continuous anisotropy field, n(r), it can be written as 
'{dn^/dri), where index /x refers to the spin 
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polarization and e^'^^ is the unit antisymmetric tensor. 
The spin current is responsible for the transmission of 
angular momentum and produces the torque rotating the 
anisotropy axis at each site to orient all of them in the 
same direction. 

In the case of magnetic rings, the first term in (14) 
can be nonvanishing, which is related to the Berry phase 
acquired by a moment moving along the ringi^. As we as- 
sume the magnetic ordering uniform, the Berry phase has 
a non-geometric origin but is determined by the topology 
of the ring^*' . In its turn, the nonvanishing Berry phase 
affects the inversion symmetry with respect to the motion 
along the ring in opposite directions. 

To specify this suggestions, let us consider now the 
mesoscopic magnetic ring in form of a thin-wall cylinder, 
with the magnetization along the z axis of cylinder, like 
described in Ref4. The hard axis is directed perpendic- 
ular to the cylinder surface. Correspondingly, it changes 
direction along a contour around the ring. We can re- 
strict ourselves by considering only the i^-component of 
momentum q corresponding to the around-ring motion 
^ (see Fig. 2) because the gauge potential in this case has 
only one nonvanishing component A^p = 1/-R, where R is 
the radius of the ring. Note that A^p is a constant, and 
it fully justifies using the Fourier transformation lead- 
ing to Eq. (8). The periodic condition imposes the re- 
quirement of momentum quantization, which reads as 

i^^pn + R = n, where n — 0, ±1, .... Then the first 

term in Eq. (14) gives us the spin current in the ring 
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where uj° = 2S [{an'^ / R^) {an'^ / R^ + Xi)Y^'^ . In the 
limit of /3 = XiR^/a ^ 1, which corresponds to the weak 
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gauge field regime, the main contribution is related to 
n ~ ^/]3 ^ 1, and we can substitute the sum in (17) by 
the integral. In this approximation, we obtain 



where 
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0.754. (19) 



As follows from (18), the spin current is zero if the 
anisotropy constant Ai — 0. By evaluating the second 
contribution of Eq. (14) into the spin current, we find 
that in this case it is small with respect to (18) in a pa- 
rameter ao/{R\/P) <^ 1. 

Using (18) we can estimate the magnitude of the me- 
chanical stress in the ring. We take Ai = 3.4 x 10"^"* erg, 
corresponding to Ai/aj] = AttM"^, M = 10^ G, oq = 
3 X 10"® cm; a/a^ 4 x IQ-^'' erg, 2R = 10 nm. Then 
we get the angular momentum X° ~ 3 x 10^^^ erg. Vary- 
ing the length of ring L by SL, we get a variation \6T'^\ = 
{X°- /L) SL. Correspondingly the force acting in the cross 
section of the ring is given by / = 6T"'/SL = and 
the stress a = f /So, where So is the cross section of the 
ring. Taking S'o = 5 x 10"^'' cm^ (5x1 nm^) we find 
a ~ 0.05 N/cm^. 

In conclusion, we demonstrated the presence of non- 
vanishing equilibrium spin currents in magnetic struc- 
tures with the homogeneous magnetization and varied 
anisotropy axis. One of the most important consequences 
of this effect is related to the correct definition of the 
non-equilibrium spin current and interpretation of spin 



current measurementsi^iiiii^. The existing attempts to 
exclude the equilibrium spin current have been concen- 
trated only to a part responsible for the torque acting 
on spins. In the coUinear magnetic system, this part of 
torque is zero. Nevertheless, as we demonstrated, the 
spin current can also transmit the angular momentum 
rotating the anisotropy axes like in the case of mag- 
netic rings. Similar to the magneto-electric effect in mag- 
netically inhomogeneous systems^, the mechanical stress 
in anisotropy-inhomogeneous systems can be called the 
magneto-mechanical effect. 

Our estimation of the stress suggests that the corre- 
sponding deformations are too small to be measured ex- 
perimentally. However, the spin current in the ring can 
be observed due to the electric polarization, as it was dis- 
cussed recently^i^iii. It should be noted that the com- 
bination of magnetoelectric and magnetoelastic effects in 
the ferromagnet can be connected with the physics of fer- 
roic materials, which are in scope of great activity nowi^. 

We calculated the spin current at T = 0, when the 
excitation of real magnons can be neglected. The ef- 
fect of temperature is twofold. Due to the excitation of 
magnons, one can expect an additional contribution to 
the spin current, related to the Berry phase from the 
Dirac stringiS. On the other hand, the magnon decoher- 
ence increases with the temperature due to the magnon- 
phonon interaction, and it suppresses the spin current. 
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